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EXISTENCE OF TRAVELLING PULSES IN A NEURAL MODEL 


S. P. HASTINGS. 

This paper is dedicated to the memory of J. Bryce McLeod, 1929-2014, a dear friend and 

inspiring collaborator. 


Abstract. In 1992 G. B. Ermentrout and J. B. McLeod published a land¬ 
mark study of travelling wavefronts for a differential-integral equation model 
of a neural network. Since then a number of authors have extended the model 
by adding an additional equation for a “recovery variable”, thus allowing the 
possibility of travelling pulse type solutions. In a recent paper G. Faye gave 
perhaps the first rigorous proof of the existence (and stability) of a travelling 
pulse solution for a model of this type, treating a simplified version of equa¬ 
tions originally developed by Kilpatrick and Bressloff. The excitatory weight 
function J used in this work allowed the system to be reduced to a set of 
four coupled ODEs, and a specific firing rate function S, with parameters, was 
considered. The method of geometric singular perturbation was employed, 
together with blow-ups. In this paper we extend Faye’s results on existence 
by dropping one of his key hypotheses, proving the existence of pulses at at 
least two different speeds, and in a sense, allowing a wider range of the small 
parameter in the problem. The proofs are classical, and self-contained aside 
from standard ode material. 


1. Introduction 


In this paper we consider the following system of equations: 

= -u (x, t) + J{x-y)q (y, t) S (u (y, t)) dy 
i = 1 - q 0, t) - pq (x, t) S ( u (x, t)) 
where J is a normalized exponential 


( 1 . 1 ) 


( 1 . 2 ) 


J(x)=° 2 e-^ 


and the “firing rate” function S is given by 


(1.3) 


1 


S (it) = -- r , -r 

V ’ 1 + e A(K-u) 


for certain positive parameters e, A, b, k. and j3. The variable u is the synaptic 
input current for a neural network with synaptic depression, the effect of which is 
represented by the scaling factor q. These equations were proposed and studied by 
G. Faye in [3]. 

The Faye model is a simplified version of one first introduced by Kilpatrick 
and Bressloff in M- These authors included a variable and equation to allow for 
spike frequency adaptation. However they show by numerical computation that 
adaptation has little effect on the resulting waves. Faye dropped the adaptation 
equation and variable in m to get his system <& See [J and na for further 
information on the physical background of (T3)- 
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In 3] the author proves two interesting results about the system B namely 
the existence of a “travelling pulse” solution and the stability of this solution. A 
travelling pulse solution of (1.1) is a non-constant solution (u,q) of the form 

(u (x + ct) ,q(x + ct)) 


such that both lim^oo (u (s), q (s)) and lim s _ ) ._ 00 (it (s), q (s)) exist and these lim¬ 
its are equal. In this paper we are interested in the existence of values of c for which 
) has such a solution. As we describe briefly below, using leads to a set 
of four ode’s in which c is a parameter. To show that a travelling pulse exists for 
some c > 0, Faye uses the theory of geometric singular perturbation initiated by 
Fenicliel in (Bj and extended by Jones and Kopell in [13]. The “blowup” method is 
also employed [2]. 

Here we extend the existence result in [4] in several ways. We show that for 
sufficiently small e > 0 there are at least two travelling pulses, hence a “fast” 
pulse and a “slow” pulse, for speeds c* > c* > 0. Also, we remove an important 
hypothesis used in [3], one which can only be verified by numerical integration of 
a related ode system. (This hypothesis is stated and discussed in Section]!]) Our 
proof is for a general class of firing functions which includes the specific S for which 
Faye states his theorem. 

Further, we use a method which allows, in some sense, a larger range of e than 
seems possible with geometric perturbation. This will be made precise in the state¬ 
ments of our theorems. We believe, based on our past experience with a similar 
problem, that it is feasible to check existence rigorously for particular positive val¬ 
ues of £ > 0, using precise numerical analysis based on interval arithmetic, but we 
have not carried out such a check. This will be explained further in Section [4] 

We now mention two well-known predecessors of the Kilpatrick-Bressloff and 
Faye models. In 1992, Ermentrout and McLeod studied the equation 

/»oo 

(1.4) Ut = -u+ J (x - y) S (u(t,y))dy. 


( 1 . 


As above, S is positive, bounded, and increasing. Since there is no feedback 
in the equation, (1.41 supports only traveling fronts, where u is monotone. In 
the landmark paper [3] Ermentrout and McLeod proved the existence of fronts 
for a wide variety of symmetric positive weight functions J and firing rates S. 
(Their work applied to a more general equation ) Subsequently, in (TB], Pinto and 
Ermentrout introduced the needed negative feedback in order to get pulses. Their 
system is 


(1.5) 


du{x,t) _ 

dt 


= -u-v + J (x -y)S(u (f, y)) dy 

1 dq(x,t ) 


They analyzed this system primarily for the case S ( u ) = H (u — k) where H is 
the Heaviside function and n is a constant representing a firing threshold. While 
some partial results have been obtained recently by Sclieel and Faye (see Section 
[!]), we are not aware of any existence proof for pulses which covers all reasonable 
smooth functions S. We discuss what we mean by “reasonable” in Section [4] where 
we also indicate why our method does not appear to apply to this model, and why 
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we expe ct that (1.5) supports a richer family of bounded traveling waves than exist 

for ([ITT]). [] 


2. Statement of results 


Travelling pulse solutions of (1.1) with (1.2) are shown to satisfy a system of 
ode’s by letting v ( s ) = f^° \e~ b ^ s ~ T ^q (t) J (it (r)) dr and computing w = v' and 
w/0 We find that 


( 2 . 1 ) 


v' = w 

w' = b * 2 3 (v — qS (it)) 


q' = e -(l-q-/3qS(u)). 

We will denote solutions of this system by p = (it, v, w, q) , and we look for values 
of c for which there is a non-constant solution p such that p (oo) and p (—oo) both 


exist and are equal. The orbit of such a solution of (2.1) is called “homoclinic”. 
In the language of dynamical systems, (u (x + ct) ,q(x + ct )) is a pulse solution of 


(1.1) if and only if the orbit of p is homoclinic. 

We make the following assumptions on S. 

Condition 1. The function S is positive, increasing, bounded, and has a continuous 
first derivative S'. 

Condition 2. The function h (it) = has one local maximum followed by one 
local minimum, and no other critical points. 


Condition 3. S is such that the system (2.1) has exactly one equilibrium point, 
say p 0 = (u o ,u o ,0,q 0 ). 

Condition 4. The function S is also such that the “fast” system 


it = 


(2.2) v' — w 

w' = b 2 (v — q 0 S (it)) 

has three equilibrium points, (ito,ito,0), (u m ,it m ,0), and (u +, u+ , 0), with uq < 
Um < u+- 

Condition 5. 

r u + 

/ (qoS (it) — it) dii > 0. 

du 0 

For convenience we will assume that 0 < S' < 1 on (— 00 , 00 ). Then Conditions 
l]|4 imply that 0 < qo < 1, uq > 0, and it+ < 1. 

We will denote solutions of (2.2) by r = (it, v, w). The local minimum of h will be 
denoted by Uk nee . In [3j specific ranges of k and A are given so that these conditions 
are satisfied by the function given in (1.3). In Figure |Tj w e sh ow the graphs of h, 
j-jgg (the q nullcline), and q = qo, when S is given by (1.3). We use the same 
parameter values as were chosen for illustration in 
We can now state our first main result. 


J We thank the anonymous referee, who helped improve the presentation and who asked several 
questions which should have been answered in the original version. 

2 It is not necessary to discuss Fourier transforms, as is usually done here. 

3 A = 20, teas 0.22 ,0 = 5, b = 4.5 
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Figure 1. graphs of h, 1+ 1 3 g , and q = go 


Theorem 1. If Conditions 0-0 are satisfied, and e is positive 
small, then there are at least two positive values of c, say c* > c*, 
has a non-constant solution p satisfying 


and sufficiently 
such that (2.1) 


lim p(t) = lim p(t) = po- 

t—f—oo t—foo 


In order to state our remaining theorems it is convenient first to give some basic 


information about the fast system, (2.2). We state this information as a pair of 


lemmas, which will be used in proving our Theorems. Their proofs are given in the 
Appendix. 

Lemma 1. If Conditions [3] and [3] are satisfied, then for each c > 0 the equilibrium 
point (u o, uq, 0 ) of (2.2fy is a saddle point, with a one dimensional unstable manifold 
Uq iC and a two dimensional stable manifold Sq, c - There is, for each c > 0 , a unique 
solution r c = (no,cWo,o R’o.c) of (2.2) with r c ( t ) £ Uo tC for all t and satisfying the 
conditions 


(2.3) 


U0,C (0) = Um 
wq :C > 0 on (—oo, 0 ]. 

Further, there is a unique c = Cq > 0 such that wo.cj > 0 on R and 

lim r c * (t) = (u+,u+, 0 ). 

£->-oo u 


In other words, the branch U(f c „ of llo lC * pointing into the positive octant u > 
uq, v > uq, w > 0 is a heteroclinic orbit connecting (uq, uq, 0) to (u+, u+, 0) . Also, 
»o,c 5 > 0 on (—oo, oo), which implies that Uq c . > 0 and u' 0 c » > 0. This solution 
is called a “front” for ( |2.2[ ). A front for ( |2.2[ ) can be characterized as a solution of 
this equation which exists on (—oo, oo), is nonconstant and bounded, and satisfies 
u' c > 0 on (—oo, oo). 


Lemma 2. If c > Cq then iCo, c > 0 on R, and 

lim Mo c (t) = lim Vq c (t) = oo. 

t—f oo ’ t—f OO 

If 0 < c < Cq then u>o, c is initially positive and has a unique zero. Also, u' 0 c has a 
unique zero, and 

lim uq c ( t) = lim uo, c (t) = —oo. 

t—f oo ’ t—f OO ’ 

If c £ (0, c(5) and t\ (c) is the zero ofu' 0 c , where uq jC is a maximum, then u'f c {t\ (c)) < 
0 and lim c _^ c »- r c [ti (c)) = (u+,u+,0). 
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Suppose finally that for some C\ £ (0, Cg) , wq, Ci > 0 on an interval (—oo, T (ci)) , 
w 0 ,ci (T (ci)) = 0 , and 

Vo,a (T (ci)) > q 0 S ( Uknee) ■ 


Then for any c £ [c\, Cg), 

uo,c (h (c)) > u k nee • 

Remark 1 . We conjecture that the condition 

vo,a (t (ci)) > u knee 

would imply the same conclusion, but we have not been able to prove this. 


The positive number Cg defined in Lemma [l] plays an important role throughout 
this paper. 

Theorem 2 . Suppose that Conditions [7]-[3| are satisfied. Suppose also that there is 
a ci £ (0,Cg), such that ifT (ci) is the unique zero ofwo, Cl (which exists by Lemma 
[ 1 |), then 

(2.4) v 0 , Cl ( T (ci)) > q 0 S ( u kne e ) • 

Assume as well that for some e > 0 there is a solution p e , Cl = (u e , Cl , u e , Cl ,w e , Cl , g e , Cl 
of j[2.1\ ) with c = c i which has the following properties: 

(i) , lim p e , Cl ( t ) = po 

t—y —oo 


.... J u' e C1 > 0 on some interval (—oo,ti) and u” Cl (ti) < 0 

^ \ u' e Ci < 0 on some interval (ti, ^ 3 ] and u E , Cl (t$) = 0 

Then for the given e there are two values of c, say c* £ (0, Ci) and c* £ (ci,Cg) 
such that ( 2.1\ ) has a homoclinic orbit. 

Figure [2] below includes a graph of the orbit of a solution satisfying ([ 1 ]) and © 
projected onto the (u, q) plane, with the points ti and 1 3 marked (as well as an 
additional point £2 which is explained later). The other solution shown in that 
figure satisfies ([ 1 ]) but not ([TI]) . 

Theorem |T| is implied by Theorem [2] and the following result. 

Theorem 3. If Conditions [7] [3] are satisfied then there is a ci satisfying the condi¬ 
tions in the second sentence of Theorem [1| Further, with this ci, if £ is sufficiently 
small, then the solution p e , Cl of \2. /[ ) with c = c\ satisfies Q) and Q of Theorem 

m 


It will follow from the proofs of these results that as e —> 0, c* —> Cg. The 
following result is all we have proved about the asymptotic behavior of c*. 


Theorem 4. 


lim c* = 0 . 

£—>-0 

However there is an M > 0 independent of e such that if there is a homoclinic orbit 
for c = c* then f- < M. 
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Remark 2. For a given pair (e, C\), the hypotheses of Theor em |D| can he verified 
by checking one solution of (2(B) at c = C\ and one solution of (2.1 ), with the given 


£ and c = Ci. For the specific model considered in H], standard numerical analy¬ 
sis (non-rigorous) easily finds specific values of (e,ci) where these hypotheses are 
apparently safis/ieiQ In the discussion section we describe how this could, in prin¬ 
ciple, be checked rigorously using uniform asymptotic analysis near the equilibrium 
points of f 2.2\ ) and \2.1\ ), and then a rigorous numerical ode solver (using interval 
arithmetic) over two compact intervals. We cite a papejQ where a similar procedure 
was followed successfully, but we have not attempted it here. 

Remark 3. In [3] only one homoclinic solution is found, and there is an extra 
hypothesis about the system (2.2). (Hypotheses 3.1) As far as we know, this hy¬ 
pothesis can only be checked by numerically solving the system (2.2) We discuss 
this further in Section^ 

3. Proof of Theorem [2] 


3.1. The fast pulse. We need two simple preliminary results about the behavior 
of solutions. 


Proposition 1. For any e > 0 the regions jv < 0,u) < 0, < q < lj 


and 


ly > 1 , w > 0 , yyTg < q < 1 j are positively invariant open sets for the system 


Proof. We are assuming that 0 < S (it) < 1 for all u. Hence, q' > 0 if q < 
and q' < 0 if q > 1. Therefore jyyrg < q < 1 j is positively invariant. Further, if 
yqyg < q < 1 then v" = w' < 0 if v < 0 and w' > 0 if v > 1. The result follows. □ 

Note as well that because S is bounded, all solutions of (2.11 exist on R = 
(—oo, oo). 


Proposition 2. Ifp = (it, v, w, q) is a solution of (2.1), and u (t) > Uknee for some 
t, then either q' (t) <0 or q(t) < h (uknee) ■ 


Proof. This follows from Condition [3j which implies that the graph of the decreasing 
function q = 1+l ^ s ^ in the (q, it) plane, where q' = 0 , passes under the point 
(^fcneet ^ (See Figure^) EH 

In the first, and longest, part of the proof of Theorem [2] we show that there is a 
“fast” pulse, with speed c*(e) which tends to Cq as e tends to zero. In the second 
part we look for a “slow” pulse, with a speed c*(e) which tends to zero as e tends 
to zero. 

We will show that for any possible homoclinic orbit, u > 0. We look for ho- 
moclinic orbits such that, as well, q < qo in (— 00 , 00 ). In searching for the 
fast solution we will consider for each c > 0 a certain uniquely defined solution 


'For the parameter values used by Faye, a standard ode solver suggests that (e, ci) = (.005, .34) 

satisfies the conditions in Theorem [2] If the conjecture in Remark^ is true then it appears that 
(s, ci) = (.05, .2) would work. 

5 (on the Lorenze equations) 

6 It apears to us that because of the degeneracy at the knee, it would be harder to verify this 
hypothesis rigorously than to do the same for (i) and (ii). 
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Pc = ( u c , v c , w c , q c ) such that p c {— oo) = po ■ We will show that there is a nonempty 
bounded set of positive values of c, called A(e), such that, among other proper¬ 
ties of Pc, either q c exceeds go at some point, or u c becomes negative. We then 
examine the behavior of p c *( E ) where c*(e) = supA(e). The goal is to show that 
Pc*(e){°°) = Po- This is done be eliminating all the other possible behaviors of 
p c * ( £ ), often by showing that a particular behavior implies that all values of c close 
to c*(e) are not in A(e). 


The following result is basic to our analysis of the full system (2.1 1 . The proof 
is routine and again left to the appendix. 


Lemma 3. Suppose that Conditionshold, and let po = (ito, ito, 0, go) be the 
unique equilibrium point of (2.1). Then for any e > 0 and c > 0 the system {2. A) 
has a one dimensional unstable manifold atpo, sayU e C i with branch U(~ c starting in 
the region {u > uq,v > ito, w > 0 , } . Ifp E , c = (it e , c > v e>c , w e>c , q e , c ) Is a solution lying 
on this manifold, then for large negative t , uq < u e>c (t) < v e , c (i) and w e>c (t) > 0. 
Also, go, c = Qo, while if £ > 0 then q E}C (t) < go for large negative t. The invariant 
manifold Uf c depends continuously on (e,c) in e > 0, c > 0. (The meaning of 
continuity here is made clear in the text below.) Finally, if Ai (c, e) is the positive 
eigenvalue of the linearization of (2.1) around po, then Ai (c,e) > Ai (c, 0 ) for each 
c > 0 and e > 0 . 


The following proposition follows trivially from (2.1l and will be used a number 
of times, often without specific mention. 


Proposition 3. 

v f 

If v! = 0 then u" = —. 

c 

v” 

If u' = it" = 0 then v!" = —. 

c 

b 2 

If u' = u" = u'" = 0 then u lv = —-q S (it). 

If q' = 0 then q" = — f3qS' (u) u' 
c 

If w' = 0 then v'" = w" = b 2 {v' — q'S (it) — qS' (it) it'). 

If q' = u' = 0 then q = 0 and q" = — j3qS' (it) u" = — ^PqS 1 (it) v . 

c c 2 

If q = v! = v' = 0 then q lv = —tt f3qS r (u)w r = —^ flqS' (it) v" = — j3qS' (it) u". 

c A c 2 c 

We use the fourth item in this list to prove 

Lemma 4. For any e > 0 and c > 0, if p is a solution on Uf c and u' > 0 on an 
interval (— oo,t], then q' < 0 on (—oo, r). 

Proof. If it' never changes sign, let a denote oo. Otherwise, suppose that it' first 
changes sign at a. If q' (r) = 0 for some r < a and r is the first zero of q', 
then q" (t) > 0, and by the fourth item of Proposition [3j it' (r) < 0. From the 
definitions of a and r, v! (r) = 0 and so q" (t) = 0. Since u' does not change sign 
at r, it" (t) = 0 and so q'" (r) = 0. Hence at r, 

u' = it" = q' = q" = q'" = 0 . 
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If u'"(t) = 0 then p(r) is an equilibrium point, a contradiction. If v!" (r) < 0 
then r is a local maximum of u !, which is inconsistent with the assumption that 
v! > 0 on (—oo,(j]. Hence u'" (r) > 0. But then q lv (r) < 0. This again implies 
that q' > 0 on some interval to the left of r, contradicting the definition of r. This 
completes the proof of Lemma [4] □ 

Lemma 5. If p = ( u,v,w,q) is a solution on U+ c then w > 0 on an interval 
(—oo, t] with u (t) = u m . 

Proof. Observe that h ( u ) > go for uq < u < u m . It follows that if uq < u < u m 
and g < go on an interval (—oo,f), then w' > 0 on this interval. Hence w > 0 as 
long as uq < u < u m and q < go- (That is, if uq < u < u m and g < go on (—oo, f], 
then w > 0 on this interval.) Since u' > 0 as long as w = v' > 0, Lemma [4] implies 
that w > 0 as long as Uq < u < u m , proving Lemma [5j □ 


Hence the conditions u 
solution 


( 0 ) = u m and w > 0 on (—oo, 
Pe,c = Qe,c) 


0 ] 


determine a unique 


on U+ c . 
Let 


Cl = < (u,v,w,q) | 0 < u < 1 , 0 < v < 1 , 


1 + /? 


< q < 1 


Since (uo,Mo,0, go) £ H, it follows from Proposition [T] that if Uf~ c is a homoclinic 
orbit, then it lies entirely in II. 

In the rest of this subsection let ci and e be chosen as in Theorem [2] 


Lemma 6. If c > ci, then either w EjC > 0 on R, or at the first zero of w EjC , 
n El c P I^/cnee* 


Proof. By the hypotheses on Ci in Theorem [2} Lemma [2] implies that 

uo,c (T (c)) > u kn ee. 

To extend this to e > 0 a comparison result is needed. Let p = (it, v, w , g) = p EtC . 
Lemma [4] implies that if £ > 0 then g < go on any interval (—oo,t] where w > 0, 
since in such an interval v' > 0 and v! > 0. Also, as long as w > 0 we can consider 
rt, w, and q as functions of v. Say that u = U ( v ), w = W (n), and q = Q(v). Then 

TT'(v) - v ~ u ^ v ) 

/"oil U [ v > - cW(v) 

_ P(v-Q(v)S(U(y)) 
vv ( v ) ~ W(v) 

We compare w = W ( v ) with the solution when e = 0. Let p\ = po, c - Then we 
can write u\ = U\ (fi), w\ = W\ (hi), and q — go- The equations become 


(3.2) 


U[{v) 
W{(v) = 


_ v—Ui(v) 

cW\{v) 

b 2 {v-q 0 S{Ui{v)) 

IVi(») 


Since Ai (c, e) > Ai ( c,0) (Lemma [3]), it is seen by considering eigenvectors of the 
linearization of ( 2 . 1 ) around pcQ that for v sufficiently close to uq (i.e. for large 


7 The relevant matrix B is given in Appendix B. 
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\ W (v) > Wi ( V) 


If, at some first v, one of these inequalities should fail while the other still holds, 
then a contradiction results from comparing © and ( |3.2[ ), because q < q 0 and S 
is increasing. For example, if U(i>) = U\(v) and W(v) > W\(v) > 0, then ( |3.1| ) and 
(3.2) imply that U'(v) < U[(v), a contradiction because U < U\ on (uq,v). Also, 
if W(v) = W\(v) > 0 and U(v) < U\(v) then W’(v) > W[(v), since q < qo as long 
as w > 0. This is also a contradiction of the definition of v. 

If both inequalities fail at the same v, then there is still a contradiction because 
q < q 0 . Hence, if W\(v) > 0 for u 0 < v < v then (|3.3| holds in this interval. 


This implies that for any c € (ci,Cg), if w E}C has a first zero at T (e,c), then 
u £iC (T (e, c)) > vo, c (T (0, c)). In the proof of Lemma[2]it is shown that Vq :C (T (0, c)) > 
vo,a (T (0, ci)), and combining these shows that if u e , c (T (e, c)) < Uknee and vq, Ci (T (0, Ci)) > 
qoS (uknee) then w’ e c (T (e, c)) > 0. This contradiction completes the proof of 
Lemma [ 6 ] □ 


Lemma 7. If c > Cq, then v e , c > 0, w e , c > 0, u' e c > 0, and u e , c —> oo as t —» oo. 
Proof. This follows from Lemma [2] and the comparison used to prove Lemma [ 6 ] □ 

We are now ready to apply a “shooting” argument to obtain the fast pulse. 
Still with £ as in Theorem [2j for each c > 0 let p = p e , c and set 


A = {c > Ci | There exist t\, t 2 , and to such that 0 < f | < t -2 < to and 

if p = p e , c then v! > 0 on [0, ti), v! (ti) = 0, u (< 2 ) = uq, and either u (to) = 0 or q (to) = qo- 
Further, u" (ti) <0, u' < 0 on (ti, < 2 ] and u < uq on (t. 2 ,to}.} 


(See Figures [2]) 




Figure 2. Each figure shows one solution with c £ A and one 
solution with c ^ A. 


Lemma 8 . A is an open subset of the half line c> C\. 
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Proof. Suppose that c £ A, and choose £3 = £3 (c) as in the definition of A. Note 
from ( 2.1 ) that if u c (£3) = 0 then there is a r < £3 such that v (t) = 0, v' (t) < 0. 
Also, v" < 0 if v < 0. Hence v (£3) < 0 and u' (£3) < 0. 

Also, ( 2 . 1 ) implies that if ^(£ 3 ) = go and u{tz) < uo then g'(£ 3 ) > 0. Since p c (t ) 
is a smooth function of c, uniformly for t in, say, (— 00 , £3 (c*) + 1 ], it follows that for 
c in some neighborhood of c*, ti (c) is defined for i = 1, ..3 and all the inequalities 
in the definition of A continue to hold, so that this neighborhood lies in A. This 
proves Lemma [Sj □ 


The hypotheses of Theorem [2] imply that C\ £ A, while by Lemma [7J if c > Cq, 
then c (j A. The numbers £j depend on c, and when we need to emphasize this we 
will denote them by £j (c), for i = 1 , 2 ,3. 

We now let c* = sup A. (This is finite, by Lemma]?]) 

Lemma 9. With e as in Theorem Q is a homoclinic orbit of 

Proof. The proof depends on the fact that c*is a boundary point of A and lies 
in (ci,Cq) . In Figure 4 we show several graphs of (w, q) which, if they occurred 
for p SiC , would suggest (without quite implying) that c was on the boundary of 
A. We must eliminate these and some other possibilities, and this will imply that 
p EtC * ( 00 ) = (w 0 , wo, 0, go) ■ The reader may want to review the definition of A before 
examining these figures. 


' 2 . 1 ). 




Figure 3. Several possibilities if c is on the boundary of A 


We need additional lemmas. Recalling that e was chosen in the statement of [2j 
we will now drop the e-dependence of p E}C and its components from our notation, 
writing p c and A. When the dependence of p c on c is not crucial to an argument 
we will use u, v, w, q for its components. 


Lemma 10. Suppose that p = (w, v, w, q ) is a non-constant solution of (2.1 ) satis¬ 
fying one of the following sets of conditions at some r : 


(i) u (r) = 0 , v' (r) = 0 , 

(ii) u' (r) = 0, v' (r) < 0, 

(iii) v! (t) = 0 , v' (r) < 0 , 

(iv) u' (r) = 0, v' (t) < 0 , 


w' (r) < 0 , q (t) < 0 , u (t) > wo 
w' (r) = 0 , q' (t) = 0 , u (r) > w 0 
q' (t) > 0 , w (r) = w 0 . 
w' (t) > 0 , q' (r) > 0 , w (t) > u 0 


Then p (r) ^W+ c ,. 
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Proof. Suppose that (i) holds. Then w' (t) and q' (r) cannot both vanish. If 
w' (t) = 0 and q' (r) < 0 then w" (r) = — b 2 q' (r) S' (u (r)) > 0. Hence in some 
interval (r — 5, t), 


(3.4) w' < 0 and q < 0. 

If q' (r) = 0 and w' (r) < 0 then 


q” (t) = --^q(r) S' {u{t)) u' (t) =0 

q'" (t) = -\()q (t) S' (u (t)) v' (r) = 0 

q (zv) (t) = ~\pq (t) S' (u (t)) w' (t) > 0 . 
& 


Once again we see that (3.4) holds on some interval 
Consider the “backward” system satisfied by P ( s ) 
then 


(t - S,t). 

= P (t-s).1£P=(U,V,W,Q) 


(3.5) 


U' = ^ 

V = -W 

W’ = b 2 (Q5 (U) - V) 
Q' = | (Q + PQS (C/)-i) 


Also, 


(3.6) 


U' (0) = 0,1/' (0) = 0, W' (0) > 0, and Q' (0) > 0. 


From (3.4) and ( |3.6[ ) it follows that on some interval 0 < s < S, 
(3.7) U' > 0, V' < 0, W' > 0 and Q' > 0. 


We claim that these inequalities hold for all s > 0. If, on the contrary, one of them 
fails at a first so > 0, then 


(3.8) U (s 0 ) >U(0),V (s 0 ) < V (0), W (s 0 ) > W (0), and Q (s 0 ) > Q (0). 


But (3.8), (3.5|, and (3.6) imply that at s 0 , all of the inequalities in (|3.7|) still hold, 


because S' > 0. This contradiction implies that U, W, and Q continue to increase, 
and V continues decrease on 0 < s < oo, and in particular, U does not tend to uq 
as s —> oo. Thus, p (r) ^ U + c ,. 

The proofs in cases (ii), (iii) and (iv) are similar and left to the reader. This 
completes the proof of Lemma [To] □ 


We now begin our study of the properties of p c *. Lemma [TO] will assist us in 
proving the following result. 


Lemma 11. The number t\(c*) is still defined, as the first zero of u' c », and 
u", (ti (c*)) < 0. Either U+ c » is homoclinic or t -2 (c*) is still defined, as the first 
zero of u — Uq- Also, ifU+ c * is not homoclinic then u' c » < 0 on (H,^]- 

Remark 4. This lemma eliminates the graph in Figure\3A\a. 


Proof. Suppose that t\ (c*) is not defined. Then u' c . > 0 on (— 00 , 00 ). Since po is 
the only equilibrium point of (2.1), this implies that for some r, y c *( T ) > u c * (r) > 1 
and v ' c , (r) > 0. Then these inequalities hold at r for nearby c, and by Proposition 
[lj v c (t) > 1 for t > t. Hence u c > 1 on [r, 00 ) and so c ^ A, contradicting the 
definition of c*. Therefore t\ (c*) is defined. 
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We now show that u", (ti) < 0. Again assume that p = p c *, and suppose that 
u" (ti) = 0. If v!" (ti) < 0 then t\ is a local maximum of u', which is not possible 
because t\ is the first zero of u’. Hence at t\, v!" = ^ (v — qS (it)) > 0. If v!" (t±) = 
0 , then at t\, 

ir" h 2 

u {iv) = _ = - b - q 'S(u) > 0, 
c c 

by Lemma [4j This implies that v! changes sign from negative to positive at ti, 
again a contradiction of the definition of t\. Hence at t±, 

u '" = h ^(. v - ( l S ( u )) > 

c* 

or q(ti) < = s™u(tl)) i s i nce u ' (£l) = 0- Also, u" > 0 in some interval 

+5). However u is bounded by 1 and does not tend to a limit above Uq. 
Therefore v! changes sign at some t > t\ (c*). Since u! > 0 on (—oo,t], v > u on 
this interval. At r, u" < 0, and so there is a point a in (ti,r) such that u" = 0 
and u'" = v ~ u = ^- < 0. Hence at a, q > But by Lemma Jl7| 

u(ti) > Uknee , and h(u) = g^jy is increasing in (ufc nee ,, oo), so q(a) > q(t 
again a contradiction of Lemma |4j We have therefore proved the first sentence of 
Lemma HH 


Lemma 12. If p = p c * t then v! < 0 as long after t\ as q' < 0. 


Proof. Suppose instead that there is a first r > t\ such that q' < 0 on (—oo, r] but 
v! (t) = 0. Then u" (r) > 0. First consider the case q' < 0 on (—oo,r]. 

If u" (t) > 0 then u' c » changes from negative to positive before q' = 0, and this 
will be true as well for c close to c*, contradicting the definition of c*. Hence suppose 
that u" (t) = v = 0. If u'" > 0, then u' has a local minimum at r, contradicting 
the definition of t. Hence 




But now the conditions in (i) of Lemma 10 are satisfied, giving a contradiction. 

We have left to consider the case that q' (r) = u' (r) = 0. Then q" (r) = 0. If 
u" (t) > 0 then q'" (r) < 0 so q' < 0 in an interval (r, t + d). Hence in this case, v! 
changes sign (from negative to positive) before q' > 0. For c close to c* there are 
two possibilities: either u' c changes sign from negative to positive before q' c > 0 , and 
so before u = wo, or else u' c < 0 in a neighborhood of r, but in a neighborhood of, 
say, r + u' c > 0 and u > u$. (See Figure |4j) Neither of these possibilities occurs 
if c £ A, so once again, c* ^ <9A, a contradiction. This proves Lemma [12} □ 


It follows that there is a first n > t\ such that q' c , (n) = 0. Also, u' c . (n) < 0, 
and (equivalently by the fourth item of Proposition [3| q ", (ti) > 0. 

Lemma 13. q' c » >0 and u' c » < 0 os long after T\ as u c * > uq. 

Proof. This lemma eliminates the graph in Figure |3T] -b. 

Let p = p c *. Since q" (n) > 0, q' > 0 and v! < 0 on some interval (n, n +5] with 
S > 0. We claim that q' > 0 on any such half-closed interval in which u' < 0. This 
follows because, by Proposition [3j q” > 0 at any point where q' = 0 and u' < 0. 

We next show that u! < 0 on any interval (iq, t\ +d] in which q' > 0 and u > uq. 
If not, then there is a first a > t\ with v! (a) = 0, q' (a) > 0 and u > Uq on (—oo, a]. 
Then u" (a) > 0. If u" (a) > 0, then v ! c » > 0 in some interval (a, a + S). In this 
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Figure 4. (it c ,g c ) for c = c* (with cusp), and two other solutions 
with c ^ A. The graph of a decreasing function is part of the q 
nullcline. 

case, for c close enough to c*, u' c changes sign after t\ but before u c < uq or else u c 
crosses uq and back again, and such c cannot lie in A, a contradiction. 

Hence, u" (a) = 0. But then, because u(a) = v (a ), 

u "' (c) = h ^{v-qS (u)) =^{u-qS ( u )). 

In the region where q' > 0 and u > uq, q < . Hence, v!" (cr) > 0 and again 

u ' c , > 0 in an interval to the right of a but before zi c » < u o, a contradiction as 
before. 

The only other possibility contradicting Lemma [13] is that there is a first r > 
7 ~i (c*) where zt c » (r) > uq and q' c * (t) = u' c , (r) = 0. We consider two cases: (a) 
<7c* (t) < go and u c * (r) > iio, and (b) g c * (r) = go, u c * (r) = uq. First consider (a). 
In an interval (r — <5, r), g'* > 0, u ' c . < 0, and u c * > uq, and so at t, if p = p c *, 
then 

u' = 0, u" >0, q' = 0, g" = 0. 

Also, 

q r " = -ip q S ' («)u" < 0. 
c 

But «" (t) > 0 is impossible because it means that even for nearby c, u' c > 0 after 
t\ but before u = u$. Therefore at r, q'" = 0 and u" = 0. Then 

g ( “) =--qS' (u)u m . 
c 

But on the nullcline q' = 0, with u > uq : u' = 0, and u" = 0, 

(3.9) u ,n = ^ ^ (v - qS (u)) = ^ (u - qS ( u )) > 0. 

This implies that u' has a local minimum at r, whereas we know that u' < 0 in 
(<i,r). This contradicts the definition of r. 

Turning to case (b), we now have that at r, 

u' = 0, q 1 = 0, q" = 0, (u, g) = (u 0 , q 0 ). 
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Thus w' (t) = 0. If u" (t) > 0 then u! > 0 to the right of r. As before, if c is 
close to c* then either u c crosses Uq twice, or p c doesn’t reach the region u < Uq 
before v! > 0, both of which mean that c ^ A. 

If u" (t) = pp = 0 then p(r) is again an equilibrium point. The third possi¬ 
bility, u" (t) = pP < 0 implies that (ii) of Lemma 10 is satisfied, and thus again 
gives a contradiction. This completes the proof of Lemma 13 □ 


If u > uq on R then u' c , < 0 and g( > 0 on (n, oo), and Up, is homoclinic. This 
proves Lemma EH □ 

Thus, for p = p c , if Up, is not homoclinic then £2 exists with u (£ 2 ) = uo and 
v! < 0 on (£ 1 , ^ 2 ]- However, there is no £3 such that u < uq on (£ 2 ,^ 3 ] and either 
u (£ 3 ) = 0 or g(t 3 ) = go, for otherwise c* G A, and this has already been ruled out. 

Therefore if £ > £2 (c*) then 0 < u c * < uq and q c , < go, for otherwise nearby 
values of c are once again not in A. If Up, is not homoclinic and p = p c * , then there 
must be a first r > t 2 (c*) with u (r) = u 0 , v! (t) = 0 , q (r) < g 0 and u" (r) < 0. 

Suppose that this is the case and also g (r) < go- (This is pictured in Figure 
3T}c.) 

Then q' (t) > 0. If u" (t) < 0, then (iii) of Lemma (10) applies and gives a 
contradiction. Hence g(r) = go- Then at r, q' = u' = w' = 0 and u" = p < 0. 


(This is pictured in Figure 3.1 -d. But this is case (ii) of Lemma 10 and so also 
impossible. 

We have established that if Up, is not homoclinic (with u > uq on R ) then for 
large £, 0 < u c * (t) < uq and g'„ >0. This is only possible if Up, is homoclinic 
(with q c , < q 0 and v! c , > 0 for large t). This proves Lemma [9] □ 

To complete the proof of Theorem [2] we look for a second homoclinic orbit, with 
c < Ci- 


3.2. The slow pulse. Again we adapt the method in m- It is stated so as to be 
useful in the proofs of Theorems and [4j as well as Theorem 

Lemma 14. There are c > 0 and M > 0, both independent of e, such that if 
0 < c < c and - > M then the solution p c remains in the region v > u on 
(— 00 , 00 ), and u crosses u = 1. 

Proof. Since the proof uses some of the easier parts of the proof of Lemma [ 2 ] it is 
included in the appendix. □ 

From here to the end of this section the parameters e and ci remain as in the 
previous subsection. Lemma[l4|implies that if we extend A to (0, Cg), with otherwise 
the same definition as above, then inf A > 0. This suggests that inf A corresponds to 
a homoclinic orbit. The problem with this argument is that the concept of “front” 
in the sense used in the method of geometric perturbation, breaks down for small 
c. The slow homoclinic orbit is not close, even up to the first zero of u' c , to the 
front found when e = 0. More precisely, our proof of the first sentence of Lemma 


Uc > Uknee- Hence we must modify our “shooting set” on the c axis. This requires 
several steps. 


11 is no longer valid, because we cannot assert that the first zero of u' c occurs with 
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Our argument from here no longer refers to a point t\ where u' c changes sign, 
but instead considers solutions such that q' c changes sign. Let 

E = {c£ (0, Ci] | There is a iq > 0 such that q' c < 0 on (—oo, iq), 

< 7 c (n) = 0 , and u' c (iq) < 0 } . 


Our argument does not require that if c £ E then u' c has only one zero in (—oo, ti) , 
though numerically this appears to be the case. 

Ei = {c £ E | q c > 0 on any interval (ti, T) in which u c > 0 and q c < go} ■ 


Recall that e and Ci were chosen so that v! has a unique zero. As in the proof of 


13 


Lemma 

Ci £ Ei 

points of E. 
Let 


this implies that g}, has a unique zero, say r Cl , and u! (t C i ) < 0. Hence 


Also, Lemma 14 shows that there is an interval (0, C 2 ) which contains no 


C 3 = sup {c < Ci | c ^ E} . 


Lemma 15. There is a iq such that q' C3 < 0 on (— 00 , iq), g ' 3 (ti) = 0, g " 3 (iq) = 0, 
and g'" (iq) < 0 . 


Proof. If g' < 0 on R then there is a a > 0 such that u C3 (a) = 1 and u' C3 (a) > 0. 
From the continuity of p c (t) with respect to c, the same is true for u c if c is 
sufficiently close to C3. In particular, again g' < 0 on (— 00 , 00 ). But then c ^ E, 
contradicting the definition of C3. 

Therefore a first iq is defined such that q' C3 (iq) = 0. Then g " 3 (iq) > 0. 
Also, by Proposition [3j g" (iq) = -f)S' (u C3 (iq)) u' C3 (iq). If g" (iq) > 0 then 
by the implicit function theorem, iq (c) is defined for nearby c as the first zero 
of g', with g" (iq (c)) > 0 and g' < 0 on (—oo,iq(c)), contradicting the defi¬ 
nition of C3. Hence g " 3 (iq) = u' c (iq) = 0. If q” 3 (iq) > 0 then iq is a local 


minimum of q' C3 , contradicting the definition of iq. If q" 3 (iq) = 0 then g, 

-^/39c 3 (n) S' ( u C3 (rj)) w' C3 (ti) , and since q' C3 (n) = 0 and g C3 (r x ) < g 0 , u C3 
0 and g“ (r x ) < 0. This implies that g} 3 > 0 on an interval (iq — <5, iq), again a 
contradiction. Hence g"( (n) < 0, completing the proof of Lemma 15 


(n) = 

(n) > 


□ 


Thus, q' C3 < 0 in some interval ( 7 q,iq + (5). This result implies that C 3 ^ E. 
However the interval (C3, ci] C E. Lemma [l5| also implies that points in (C3, ci] near 
to C 3 are not in Ei, since the corresponding solutions on IAf~ c must have a change 
of sign of q' c from positive to negative after iq (c). Let 


c» = inf {c > C 3 j c £ Ei} . 

We claim that Uff is a homoclinic orbit. 

The proof uses techniques very similar to those above. First observe that c* > C 3 
and c* £ E. Therefore iq = iq (c*) is defined as in the definition of E. Then use the 
following result. 


Lemma 16. If c £ Ei, then u' c < 0 on any interval [iq (c), iq (c) + 5] in which 
u c > u 0 . 


Proof. If u' c = 0 at some first a > iq with u c (a) > uq, then u" (a) > 0. But in the 
region where u> Uq and q' > 0 , w' is positive, and this implies that p c crosses into 
q' < 0 , a contradiction of the definition of Ei. □ 
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Corollary 1. IfU+ Ct is not homoclinic then there is at 2 > ti such that m c „ (t 2 ) = 
Mg and u' Cf < 0 on [tl^] • Further, u Ct < Mg on (£ 2 , 00 ). 


Proof. Let p = p c> . Lemma 16 implies the existence of t 2 . Suppose there is a first 
<7 > t 2 with u (a) = uq. From the definitions of Si and c*, q' > 0 on [tr, 00 ). Since 


is decreasing and qf > 0 if q < inf 


l+/3S(u 

Hence there is a p > a with 


1+/3 S 


u cannot increase indefinitely. 


u ( p ) > M 0 , v! ( p) = 0 , m" ( p ) = 

c 

w' ( p) > 0 , q' (p) > 0 . 


w(p) 

c 


< 0 


A contradiction then results from (iv) of Lemma 10 


□ 


Now apply the technique of Lemma[ 8 j including use of P rop osit ion [3| an d Lemma 
10 to show that u Ctt > 0 and g c „ < qo on (t 2 , 00 ). In particular, Lemma|lO|is used to 
show that there is no t > t 2 (in fact, no t at all) with (m c „ (t), q Ct (t)) = (ug, qo). It 
follows that on (f 2 ,oo), q' c , > 0, and so indeed, Uf~ c * is homoclinic. This completes 
the proof of Theorem [2j 


3.3. Proofs of Theorems [l] [3] and 4l As mentioned above, Theorem [T] follows 
from Theorems [2] and [3j In Theorem 3je is not fixed. Also, Ci is any number in 
( 0 , c(j), which however is fixed at this stage for the rest of this section. 

Let I = [ci, Cq + 1]. We note that the unstable manifold Uf[ c varies continuously 
with (e, c) for e > 0 and c £ I. To be more precise, if exists on (—00, T] then in 
some neighborhood of (e, c), p £jC (f) exists for —00 < t < T and is a continuous 
function of (e, c, t ) [^] 

From the third sentence of Lemma [ 2 ] it follows that Ci can be chosen in (0, Cq) 
such that if ci < c < Cg, then vo, c {T (c)) > qoS ( Uknee ), where T (c) is the unique 
zero of wo,c- This proves the first assertion of Theorem[3j We now choose Ci in this 
way. 


Lemma 17. There is an Eq >0 and a r > 0 such that if 0 < e < £q and 
c £ I = [ci,Cq + 1 ], then u e<c ( t ) = Uknee for some first t < r, and w £iC > 0 on 
(— 00 , t]. (Hence, p EtC satisfies the the first condition of Theorem [Ip Further, Eq 
can he chosen so that p s , Cl satisfies conditions |*|j and ^ of Theorem 

Proof. From the choice of c\, Lemma [2] implies that for any c > ci, there is a 8 > 0 
such that if c £ I then v' 0 c = wq^ c > S in the interval [0, r] where u m < u < Uknee- 
From cm' = v — u it follows that for some r > 0, if ci < c < Cq then mo iC = Uknee 
before t = r. The uniform continuity of Mo jC (i) in (e, c, t ), for —00 < t < r and d £ I 
in any compact interval [0, £] with e > 0, then implies the first conclusion of the 
Lemma. The remaining assertion of Lemma [TT] follows by similar arguments. □ 

We have now proved Theorems HJ |3j and [TJ in that order. To prove Theorem g] 
apply a continuity argument similar to that just used to show that for any 6 > 0 
there is an £1 such that if S < c < Cq — S and 0 < £ < £ 1 , the pair (£, c) satisfy the 
hypotheses on (£, Ci) in Theorem [ 2 ] It follows that pulses exist for some c* € {0,6) 
and some c* € (cq — <5, Cq) . But Lemma [l4| implies that f- < M. Theorem|4lfollows. 


®See the footnote at the end of the appendix for a further discussion of this point. 
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4. Discussion 


4.1. Hypothesis 3.1 of [2]. As stated earlier, there is an additional hypothesis 
in the existence result given in [5|, namely Hypothesis 3.1 in that paper. This 
hypothesis is interesting in a broader context, and we will include some comments 
on its relation to the well-known pde model of FitzHugh and Nagumo. 

To state this hypothesis we need to introduce a basic tool in the method of 
geometric perturbation, the so-called “singular” solution. The singular solution of 
(2.1) is a continuous piecewise smooth curve in l? 4 consisting of four smooth pieces. 

The first piece is the front with speed Cq found in Lemma[lJ (Recall that “fronts” 
were defined just after the statement of this lemma.) In Figure [5] the green line 
segment is the projection of the graph of the front onto the (u, q) plane. 



Figure 5. The red curve is the graph of q = 


The second piece of the singular solution is a segment of the nullcline u = qS (it) 
(with w = 0, v = u) as shown in Figure [hj It is obtained from (2.1) by letting 
p(t) = P (et ), formally setting e = 0 in the resulting system of ode’s for P = 
(U,V,W,Q ), and solving the resulting set of one differential equation and three 
algebraic equations, one of which is U — QS (U) = 0. For more information on this 
segment, and the singular solution in general, see [4]. We don’t need to say more 
about this segment here. 



Figure 6 . The first two segments of the singular solution 


To define the third part of the singular solution (crucial in Hypothesis 3 of |4]), 
we consider the fast system (2.2), but with q 0 replaced by q, as a parameter ranging 
between g m j n and go- 
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(4.1) 


[/ = b 2 (v — qS (u)) 


For each q £ {q m in, </o] there is a unique c(q) > 0 such that (4.1) has a bounded 
non-constant solution 


r c(g) — { U c(q)l V c(q)l W c(q)) ■ 

The graph of this solution is a heteroclinic orbit connecting the left and right 
branches of q = in the (it, q ) plane. There is at least one value of q £ (<?min, <Zo) 
such that c (q) = 0. This will be true for any q such that 


ru+iq) 




(qS (u) 


u) du = 0. 


For q £ (g m in,9o) and sufficiently close to q m i n , the integral above is negative, 
and whenever this is the case, the connecting heteroclinic orbit exists for some 
c (q) > 0, but unlike the front defined earlier, u c ( g ) is decreasing, from u+ (q) to 
U- (q). Such a solution is called a “back”. If q = q m - m then there is a “back” for 
any c > lim + c (q ). All this can be proved using methods from the appendix, 
or see 0. 


The third part of the singular solution of (2.1l is a “back”, at a value q = qi £ 
[9min,9o) with speed c = c(qo). There must be at least one q± for which such a 
back exists. As far as we know there is uo proof that c is a monotone function 
of q in (<7min><7o) > so possibly there could be more than one such q\. In this case 
we can require that the jump down is at the largest possible qi supporting a back 
with speed c(qg). The singular solution is said to “jump down above the knee” if 
qi > (/min* 

Numerical computations suggest that often no such qi exists in (<y m i n , <?)• In this 
case, there is still a traveling back with speed c (qo ), but it is at q\ = ( 7 m i„. 



Figure 7. Nullclines: Red, Blue; Singular solution: Green 

The fourth part of the singular solution is a slow return along the left branch of 
the nullcline. 


Hypotheses of Faye: 
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(i) The system (2.1) has a unique equilibrium point (uo, uo, 0, qo). 

(ii) If h{u) = then the equation 

9 (u) = q 0 


has exactly three solutions, uq < u m < u+, with h' (uo) > 0, h! (u m ) < 0, and 
h! (u + ) > 0. 

(Hi) 


r u + 


(q 0 S ( u) 


u) du > 0. 


(iv) The speed of any “back” with q £ (q m in,qo) less than c(qo). 


It follows that the back of the singular solution is required to be at the knee. 
This condition can only be verified by numerical integration of the fast system. 


Theorem 5. (Faye): Under hypotheses (i), (ii), (Hi), and (iv), if £ is sufficiently 
small then the system (2.1\ ) has a homoclinic orbit for at least one positive value of 
c. 


Geometrical perturbation, based on work of Fenichel and others, is a technique 
for showing that the singular solution is close to a real solution if e is sufficiently 
small. Certain “transversality” conditions can be complicated to check, requiring 
a technique called “blow-up”. 

We can contrast the existence of a back at the knee with the well-known behavior 
of the pde model of FitzHugh and Nagumo. (See [TO] for a presentation of this 
model and a proof that it has two traveling pulses. ) The fast FitzHugh-Nagumo 
pulse can be described as having a jump up, (close to a front) during which u 
increases rapidly while w is nearly zero, followed by a slow increase in w, and then 
a jump down (close to a back), with w again nearly constant (but positive), and 
u decreasing rapidly. In this case, the back occurs before ( u , w ) reaches the knee. 
See El 



Figure 8. FitzHugh-Nagumo pulse (it, v plane) 


We have done some preliminary numerical investigation to test whether it is 
possible, in the model studied in this paper and with the particular function S in 
(1.3), to adjust the parameters A and k so that the jump down occurs before w 
reaches the knee. We have not found such a pair (A, n ), but we cannot assert that 
none exists. 
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For the FitzHugh-Nagumo model, however, it is clear that the jump down is al¬ 
ways before reaching the knee. This follows because the reaction term in equations 
is a cubic polynomial, / (it) = it (1 — it) (it — a), where 0 < a < - 5 . This function is 
symmetric around its inflection point, which leads to the ’’before the knee” behavior 
of the singular solution. 

So we searched numerically for alternative functions to use for / which, while 
still “cubic like”, permit the down jump of the singular solution to be at the knee. 
We found such a function, as illustrated in Figure [9] We are not aware of a method 
which determines analytically where the downjump occurs, either for the model of 
Faye or that of FitzFlugh-Nagumo when / is asymmetric. 



Figure 9. a ’’cubic-like” function for an alternative FitzHugh- 
Nagumo type model 

The existence proof in [TOj for fast and slow homoclinic orbits of the FitzHugh- 
Nagumo system also applies to functions / such as that pictured in Figure [9] On 
the other hand it appears that the proof by geometric perturbation in this case, 
while probably still basically valid, requires a more complicated analysis because 
the downjump of the singular solution may occur at the knee. 

4.2. The slow pulse. It appears that the fast-slow analysis used to get the singular 
solution, for any of the models we have discussed, does not apply to the slow pulse. 
Hence it appears difficult to use geometric perturbation to obtain this solution. 
In [ j!5| this solution was obtained in the FitzHugh-Nagumo case using dynamical 
systems methods, but only for a sufficiently close to \. 

The analysis in m was in the region a close to |, e and c small. Indeed, if 
a = \ then there is only one pulse and it is a standing wave (c = 0). It is proved in 
na that for 1 — a positive but small there is a smooth curve (c, e (c)) corresponding 
to pulse solutions and connecting c* to c*. The speeds c* and c*, and the maximum 
of in [c*, c*} all tend to zero asa->| . 

We suspect that a similar picture holds for the model of Faye, but it is not clear 
that our analysis is able to prove this much. The FitzHugh-Nagumo condition 
that a — \ is small would be replaced here by requiring that J^ + ( q 0 S (it) — it) du 
is small. 

4.3. Model of Pinto and Ermentrout. We said in the introduction that there 
was only a partial existence result for the model in |16j . We were referring there 
to a recent paper by Faye and Scheel [5|, which uses an interesting extension of the 
geometric perturbation method to infinite dimensional spaces to handle this kind 
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of problem. The method is powerful because it allows an extension beyond the 
sorts of kernel which reduce the problem to an ode. 

Faye and Scheel remark that their paper appears to apply to the Pinto-Ermentrout 
model. This is true, but with a limitation. A key hypothesis in their paper is that 
for the singular solution, as described above, the jump down occurs above the knee. 
In a private communication Professor Ermentrout has observed that while this is 
true for the Pinto-Ermentrout model in some parameter ranges, it is also common 
for the down-jump to occur at the knee. This is why we characterized their result 
as ’’partial” for Pinto-Ermentrout. 

Unfortunately, we have not been able to make our approach work for Pinto- 
Ermentrout. The reason may be related to an important difference between (2.1) 
and the equivalent set of ode’s obtained from (1.5). The linearization of (2.1) 
around its equilibrium point has only real eigenvalues, for any e > 0 , while the 
equivalent linearization for (1.5) has complex eigenvalues for a range of positive e. 
Thus, a homoclinic orbit would oscillate around equilibrium. A few oscillations 
could occur even for the very small values of e where the eigenvalues are real. The 
final steps in our proof above clearly do not allow such oscillations. 

In [0] we showed that there was co-existence of complex roots and a homoclinic 
orbit for FitzHugh-Nagumo, and observed that work of Evans, Fenichel and Feroe 
then implied the existence of many periodic solutions and a form of chaos. This 
leads to a conjecture that the Pinto-Ermentrout model supports a richer variety 
of bounded solutions than the model of Faye. The new solutions are probably 
unstable, however, so their physical importance is unclear. 


4.4. Stability of these solutions. A local stability result for the fast solution 
was proved by Faye. His proof depends crucially on analysis of both the front and 
back of his solution, as described above. The analysis of the back is less standard 
because of the assumption that the jump down is at the knee. Here he relies on 
previous work on similar problems. Presumably if the jump were above the knee 
(where, however his existence proof is not claimed to apply), the stability analysis 
would be easier. 


4.5. Can the hypotheses of Theorem [2] be checked rigorously for a specific 

(e, Ci)? Condition (III) says that the solution p EjCl is on the unstable manifold 
of (2.1) at the equilibrium point pq. To check condition (jii]) we must follow p E)C1 
until a point where u EiCl = 0. Our proposal for doing this is based on [7], where a 
similar procedure was followed for the well known equations of Lorenz. 

Using a standard ode solver we can arrive at a conjectured value for (e, C\). To 
begin analyzing p E}Cl numerically we would expand the solutions around pq. A high 
order expansion of U+ c results in algebraic expressions which are then evaluated 
using rigorous numerical analysis based on interval arithmetic. With this technique 
one hopes to show that U+ c enters a very small box near p 0 . For the example in |Yj 
this box had a diameter of about 10 -68 . This gives us an initial estimate accurate 
to (say) 68 significant digits. 

From there, a rigorous ode solver, as described for example in [I], would be used 
to continue p e , Cl until w EtCl = 0. Whether this can be done cannot be determined 
ahead of time. One has to run the solver. The number of guaranteed accurate 
digits decreases as the integration proceeds. We then hope that some significant 
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digits would be maintained long enough to reach u = 0. It must be checked along 
the way that maxw £)Cl > u^nee- 

Based on the great sensitivity of the Lorenz equations to the initial conditions, 
we expect that this would be easier for the Faye model than it was in [7]. We 
are not aware of any proposal to try to estimate e for the method of geometric 
perturbation. 
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Appendix A. Proofs of Lemmas 0S AND [14) 


A.l. Proofs of Lemma [l] and Lemma [2} We prove these results together. The 
linearization of (2.2) around r 0 = (mq, Uq, 0) is the system Q' = AQ with 


A = 


i 


0 


-b 2 q 0 S' (w 0 ) 


l 


0 


b 2 



The characteristic polynomial of A is 

(A.l) / (X) = X 3 + r X 2 -b 2 X--( 1 - («o) 9o ), 

c c 
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Recall that h(u) = gj^j- Condition |4| implies that the equation q 0 = h (u ) has three 


solutions, u 0 < u m < u + , and by Condition [2} h! (u 0 ) > 0. It follows that 
(A.2) q 0 S'(u 0 )<l. 

Therefore / (0) < 0. Also, /' (0) = — b 2 , and both f" and f” are positive for X > 0. 
Hence A has one real positive eigenvalue. Also, / (—4) = 'frS' {u o) qo > 0, which 
implies that A has two real negative eigenvalues. 

Further, it is easily seen that there is an eigenvector corresponding to the positive 
eigenvalue of A which points into the positive octant. If r = (u, v, w) is a solution 
lying on the branch Uq c of the unstable manifold of ( 2 . 2 ) at ro, then initially, it', 
v' = w, and w' are positive. It follows from the first two equations of (2.2) that 
v > u as long as w > 0 . Also, u > qoS ( u ) for Uq < u < u m and so w' > 0 while 
u is in (uq, u m \ . Hence there is a first to such that u (to) = u m , and we can assume 
that to = 0 . We have now proved the assertions of the first and second sentences of 
Lemma |T] 

For the third sentence of Lemmajl] and for all of Lemma[2] we need a comparison 
lemma. For each c > 0, let r c = (u c ,v c ,w c ) be the unique solution of (2.2) on 
Uq c such that w c > 0 on (—oo,0] and u c (0) = u m . Suppose that w c > 0 on a 
maximal interval (—oo,T (c)), where possibly T (c) = oo. Then in (—oo,T (c)) we 
can consider u and w as functions of v, letting u c (t) = U c (v c (t)) and w c (t ) = 
W c ( v c (t)). This defines the functions U c and W c on the interval 

I r = (un, lim v c (t)), 


= (u 0 , lim 
i—fT(c)“ 


and for v in this interval, 


(A.3) 


U' c (v) = 


- U c (v) 


W' (v) = 


b 2 (v - q 0 S ( U c (v))) 


cW c (v) cw W c (v) 

Lemma 18. If d -2 > c ?2 > 0 then Id x C Id 2 - In the interval Id lt 

U d2 < U dl 

W d2 >w dl ■ 


(A.4) 


Proof. The first sentence follows by proving (A.41 on the smaller of the two intervals. 
We first show that these inequalities hold on some initial interval uq < v < uq + 6. 
This is seen by comparing unit eigenvectors corresponding to the positive eigenval¬ 
ues Ai (di) and Ai (^ 2 ) of the linearizations of (2.2) around r 0 . Suppose that for a 
particular c the eigenvector corresponding to Ai (c) is (ni (c), 712 (c), 713 (c)). Then 

, . t 7 2 (c) 

ni (c) = FT W 

^3 (c) = Ai (c) 77 2 (c) - 


Inequalities (A.4) follow near ro if Ai (cfo) > Ai (d\). For this we turn to the char¬ 
acteristic polynomial of A, given in but now denoted by / (X, c). 

It is easier to work with F = cf , noting that c > 0. The positive eigenvalue of A 
is determined by the equation 

F (Ai (c), c) = 0 

and the condition Ai (c) > 0. Then 
dF (Xi 


dx 


dAi (c) dF 

(c),c) ^ (Ai (c),c). 
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Since F (0, c) < 0, ^ (0, c) < 0 and (A', c) > 0 for X > 0, (Ai (c), c) >0. 
Also, ^ (Ai, c) = Ai (Af — b 2 ) . It follows that dX ^ > 0 if Ai < b. But 

F (b , c) = — S' ( u 0 ) q 0 > 0, 
c 

so indeed, Ai (c) < b. 

Therefore ( |A.4 ) holds on some interval {uq,uo + 6). Suppose that the first in¬ 
equality fails at a first v £ I dl , while the second holds over (uq,v]. Then at v , 
U\ = U 2 , W 2 > W±. But then, U'^ (v) < U[ (v ), a contradiction since U 2 < U\ on 
(0, t)]. A similar argument eliminates the other possibilities, using the fact that S 
is increasing, and this completes the proof of the Lemma 18 □ 

Corollary 2. IfT (d±) = oo then T (d 2 ) = oo. 

Lemma 19. If d 2 > d\, T (df) < oo, and v dl (T (di)) > qoS (ufc ne e) j then either 
T (d 2 ) = oo or u d2 (T (d 2 )) > u kne e- 


Proof. Suppose that T (d 2 ) < oo and Ud 2 (T (d 2 )) < u k nee- Lemma 18 implies that 
Vd 2 (T (d 2 )) > v dl (T (di)) and so 

w' d2 (T (d 2 )) = b 2 ( v d2 (T (d 2 )) - q 0 S ( u d2 (T (d 2 )))) > b 2 ( v d2 (T (d 2 )) - q 0 S ( u knee )) > 0, 
a contradiction because T (d 2 ) is the first zero of Wd 2 ■ □ 

Next we must show the existence of Cg. 

Lemma 20. For sufficiently large c, T (c) = oo. 


Proof. From (A.31, in the interval (—oo,T(c)), where v c is increasing, v c > u c . 
Recall that 

u 

q i \ ^ 9o 

S ( u) 

in (uq, u m ). As long as u o < U c ( v c ) < u n 
dU c v — U c 


dW r . 


and v c < 1 , 
< 


1 


cb2S ^ U c){sm- q o) C7 l(sWI)- q o) 


,S(C/ C ) 

This implies that for large c, in the interval where uq < u c < u m , w c grows 
rapidly and so, in turn, does v c .. In particular, v c > 1 before u c = u m , and this 
implies that T (c) = oo. □ 


Now we wish to show that for small c > 0, w c < 0 before v c = 1. It is in this 
step that Condition [5] is used. 

Lemma 21. There is a ui > 0 such that for any c > 0, if \w c {t) | > w and 
0 < v c (t) < 1, then |u> c | > w for t > r and v c leaves the interval (0,1). If 
w c (t) > w then v c crosses 1 , while if w c {t) < —w then v c crosses 0. 

Proof. Let w = \/2b. Since \w'\ < b 2 , if w (r) = V^b, then for s > 0, w (r + s) > 
y/2b — b 2 s, from which follows that v must leave ( 0 , 1 ) before s = □ 

Lemma 22. If 0 < w c < w on (—00, t\ then 0 < v c — u c < cw on this interval. If 
\w c \ <w on (—00, cr], then \v c — u c \ < cw on this interval. 

Proof. With r = r c , (v — u)' = w — < w — so if v — u > cw then 

(v — u)' < 0. Also, if v — u = 0 in (— 00 , r) then(u — u)' > 0. Since v — u —> 0 + as 
t — t- — 00 , the first sentence of the lemma follows and the second is similar. □ 
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Based on this lemma, we consider, in addition to (2.2), the system 


(A.5) 


V = w 


w' = b 2 (v — qoS (v)) 


This system has equilibrium points at (uq, 0), (u m , 0), and (u+, 0), and a standard 
phase plane analysis, assuming Condition pi shows that the positive branch Uq 0 of 
unstable manifold of (A.5) at (uo,0) is homoclinic. Also we consider the system 


(A. 6 ) 


u' = b 2 (v — qoS (v — c)) 


for small c. Choose c so small that this system also has three equilibrium points, 
and a homoclinic orbit based at the left most of these. This orbit entirely encloses 


the homoclinic orbit of (A.51. 


Finally we consider the system 
(A.7) 




V = w 

w' = b 2 (v — qoS (v + c)) 


For sufficiently small c this system a lso h as a homoclinic orbit. This orbit lies 
entirely inside the homoclinic orbit of (A.5). However, the lower left branch Uq 0 


of the unstable manifold of this system crosses the homoclinic orbits of (A.5) and 


(A. 6 ), and this branch will play a role below. (See Figure 12.) 



Fig ure 1 0. homoclinic or bits o f, from inner to o uter,) A.7), 
and (A. 6 ), part of U 00 for (A.7), and an orbit of (2.2) (dotted) 


From now on, (iq,u>i) 
of the systems (A.5) 


(A. 6 ) 


(v 2 ,W 2 ), and (^ 3 , 1113 ) will denote the unique solutions 
and (A.7) respectively which he on the homoclinic 
orbits of those systems and satisfy Vi (0) = u m . In each of these cases, if (v,w) 
is homoclinic then |u;| is bounded by w. This follows from the definition of w in 
Lemma 21 the results of which also apply to (A. 6 ) and ( |A.7 ), with the same proofs. 
If |w| exceeds w then p is not bounded. 

Recall that in Lemmas [l] and [2j ro iC = (uo,c> Wo,c> Wq iC ) denoted the unique solu¬ 
tion on the unstable manifold Uq^ c such that Mq jC ( 0 ) = u m and Wq >c > 0 on (— 00 , 0 ]. 
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In the rest of this proof we will denote this solution by (u, v, w) . By Lemma 22 we 
can choose c so small that if ti is the first zero of w c , then 

(A. 8 ) 

on (— 00 , t\\. 


V — C < U < V 


Lemma 23. Condition A.8 implies that if w c > 0 on (— 00 , i], then (v c ( t) ,w c (t)) 
is in the annular region between the orbit of (vi,wi) and the orbit of (v 2 ,W 2 ) ■ 

Proof. The proof is similar to the proof of Lemma [18] As long as w c > 0, u c and 
w r can be considered functions of v r . Also, 


(A.9) 


dw c 

dv r . 


b 2 (v c - q 0 S (u c )) 


By considering the eigenvalues of the linearizations of ( |A~5| as functions of c we can 
show, using (A. 2 ), that for large negative t, (v ( t ), w (t)) lies in the claimed annular 
region. Suppose that for some first r, ( v c (r), w c (r)) lies on the upper boundary of 
this region, that is, on the homoclinic orbit of (A. 6 ), at a point where w > 0. The 
slope of this homoclinic orbit at this point is 

dw 2 _ b 2 (v 2 - gpS (v 2 - c)) _ b 2 (v c (r) - q 0 S (v c (t) - c)) 

dv 2 w 2 w c (r) 

But u c (r) > v c (t) — c (since u c < v c as long as w c > 0), and since S is increasing 
and w c (t) > 0 , it follows from (A.91 that and so the curve (v c ,w c ) 

arrives at this point from outside of the annular region, contradicting the definition 
of t. In a similar manner it is shown that ( v c ,w c ) lies above the orbit of (vi,Wi) 
as long as w c > 0. This uses the bound u c < v c as long as w c > 0. 

A similar comparison shows that if t = t\ (c) is the first point where w c (t) = 0 
then v c (ti (c)) is an increasing function of c, for 0 < c < Cg, and that for c > Cg, 
w c > 0 on R. This shows the uniqueness of Cg. To complete the proof of Lemma 
[ 2 ] we show that from the first t\ where w c (ti) = 0, the curve (v c ,w c ) lies either 
to the right or belo w the orbit ( 113 ,^ 3 ), and also below the left branch of the 
unstable manifold of (A.71, at least up to the point where w c = —w (= —y/2b ). (If 
w c = —w, then, as in Lemma [2l[ v c becomes negative, which is what we are trying 
to show. See Figure 8 .) This follows by the same sort of comparison as above, now 


comparing ( v c ,w c ) with the lower half of the unstable manifold of (A.7). This is 
possible because by Lemma 22 u c < v c + c as long as —w < w c 


< 0 . 


To prove that Cg as defined in Lemma [2] exists, we note that the set of c such 
that w (si) < 0 for some Si is open, as is the set of c such that v (s 2 ) > 1 for some 
s 2 and v > 0 on (— 00 , s 2 ] . This follows because p c is a continuous function of c. 
Lemmas [20] and [23] imply that these sets are nonempty, and their definitions and 
Proposition [l] imply that they are disjoint. Since the interval (0, 00 ) is connected, 
the existence of some positive Cg which is not in either set. Its uniqueness follows 
from the Corollary to Lemma 12. 

From the definition of Cg, w c * > 0 on R. Suppose that there is an s with w c * (s) = 
0 and w c * (t) > 0 on (— 00 , s). Then w' c » ( s ) = 0, w "» (s) > 0 and 

(s) = -b 2 q 0 S' (w c3 (s)) u' c ,. (s). 

Because v', > 0 on (— 00 , s), u', ( s) > 

c O v 7 ^0 v 7 

completes the proof of Lemma [T) 


0 , giving w'f. (s) 


< 0. This contradiction 
□ 
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To complete Lemma [2] we must prove the assertions in the third and fourth 
sentences. In the third sentence, 


w 

c 


when u' = 0, and at the first zero of u', w < 0, so u" < 0. The implicit function 
theorem and the comparison dAlt imply the limit statement. 

For the last sentence of Lemma [2j it suffices to prove that Uo jC (T (c)) > Uknee. 
Suppose instead that tio jC (T (c)) < Uknee. Since c > ci, Lemma 18 and the hy¬ 
potheses of Lemma [2] imply that 

Vo, c {T (c)) > Vo tCl (T (ci)) > q 0 S ( u knee ) 

> qoS (uo,c (T (c))). 

Hence w' (T (c)) > 0, which contradicts the definition of T (c). This completes the 
proof of Lemma [2] 

A.2. Proof of Lemma 1141 


Proof. This result is about system (£3. However the argument in Lemma [21~ 
initially about system (2.2), applies equally well to (2.1), so if w c increases mono 


tonically to above w then v crosses 1, followed by u. Hence we can assume that if 
si is the first zero, if any, of w c , then w c < w on (0,si). As earlier in obtaining 
( A.8|), it follows that if p = p c , then 


(A.10) 

on (0, Si). Therefore, 

(AH) 


V > u > v — cw 


lim [u — v) = 0 

c->0+ 


uniformly on (0, si] and for e > 0. 
Also, 


(A.12) 


lim I q — 

* — >-00 V 1 


■ps{ u) 


= 0 , 


uniformly on (—oo, si] and for c > 0, e > 0. This is proved by the same argument 
which lead to (A.ll), . 


Now consider the equation obtained from (2.1) by formally setting c = 0 in (2.1), 
namely 


(A.13) 

where 


= b 2 Z(v), 


Z (v) = v — 


S(v) 


1 + PS ( v )' 


Because (2.1) has only one equilibrium point, Z (v) > 0 if v > 0. Hence there is 
a t > 0 such that if for some t, u m < v (t) < 1 and v' ( t ) > 0, then v (t + t) > 1. 
(Here t is independent of the particular solution involved.) Lemma 
from (A.ll) and (A.12). 


14 


then follows 
□ 
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Appendix B. Proof of Lemma [3] 


Proof. Suppose that the linearization of (2.1) around p 0 is P' = BP. Then 

( 


(B.l) 


B = 


0 


-b 2 q 0 S' (u 0 ) 

V -\P0oS' (u 0 ) 

The characteristic polynomial of B is 


0 
b 2 
0 


\ 


~b 2 S (u 0 ) 
-|(1 + PS(uo)) ) 


(B.2) 


1 


g(X) = X 4 + - (1 + e (PS (uo) + 1)) X 3 + 

c 


3- (q 0 S r (m 0 ) - 1 - e (/ 3S ( u 0 ) + 1)) X 

c 

+ K? ( QoS' (u 0 ) - 1 - (3S (u 0 )) 


—b 2 H— (/3S (u 0 ) + 1) I X 


While proving Lemma [l] we showed that if £ = 0, then one of the non-zero 
eigenvalues of B is positive and two are real and negative. We also saw that 
q 0 S’ (uq) < 1, and therefore, detH < 0 if e > 0. Since the trace of B is also 
negative, if e > 0 then B has either one or three eigenvalues with negative real 
part, and for sufficiently small - it has three, all of which are real. In fact, since 
g (0) < 0, g' (0) < 0, and g'" (X) > 0 if X > 0, B has exactly one real positive 
eigenvalue for every (e,c) in the positive quadrant e > 0,c > 0. For each c > 0, as 
£ increases the other roots of g remain in the left hand plane unless, for some £, 
two of them are pure imaginary. Consideration of the characteristic polynomial in 
this case (one negative, one positive, and two pure imaginary roots) shows that the 
coefficients of X and A' 3 have the same sign. This is not the case with g , because 
the coefficient of X 3 is positive and the coefficient of X is negative. 

Hence, as asserted in Lemma |3j the unstable manifold U EiC of (2.1) at po is 


one dimensional. Further, because q 0 S r (it 0 ) < 1, it follows from (B.l) that if 
p = (/zi, p 2 , M 3 , /m) is the unit eigenvector of B with pi > 0, then p .2 > 0, and 


p -3 > 0. Also (B.l) implies that if £ = 0 then = 0 and if e > 0 then /m < 0. The 
claimed behavior for large negative t of solutions on U e ^ c follows. The continuity of 
for c a () follows from Theorem 6.1 in chapter 6 in the text of Hartman, ffU-Q 
The final assertion of the lemma, that Ai (c, e) > Ai (c, 0) if c > 0 and £ > 0 
follows by writing the characteristic polynomial of B in the form 

g (*) = Xf (X) + - (1 + (3S (u„)) / PO - £ -b 2 S (u 0 ) f3S' (u 0 ) q 0 . 

c c 

We see that since /(Ai(c, 0)) = 0, g(Ai (c, 0)) < 0. Hence Ai(c,e) > Ai(c, 0), 
completing the proof of Lemma [3) □ 


^Our terminology is different from Hartman’s because we define stable and unstable manifolds 
even when e = 0. In this case the unstable manifold lA^ c -> all we need, is the set of all solutions 
p (i) which tend to po at an exponential rate as t —»• — oo. To obtain the desired continuity of IA+ C 
with respect to e and c, apply Hartman’s theorem to \2.2\ augmented with equations s' = 0 and 
c' = 0. This is the closest we come to center manifolds in our approach. 











